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O ■ ABSTRACT 

Tj- ! We discuss the imphcations of the more than 50 microlensing events seen by the EROS, 



:z; 



MACHO, and OGLE collaborations for the composition of the halo of our galaxy. The event 
rates indicate that the halo mass fraction in MACHO's is less than 30%, consis 
expectations for a universe whose primary component is cold dark matter. We cau 
the uncertainties are such that a larger MACHO fraction cannot yet be excluded. 



In 1986 Paczynski suggested microlensing as a probe of dark (or very faint) stars in our 
galaxy |jl| (referred to generically as MAssive Compact Halo Objects, or MACHOs). In the 
past year more than 50 microlensing events have been reported. (For microlensing the two 
images are too close to be resolved; instead, the combined light leads to an achromatic, 
time-symmetric brightening.) The EROS collaboration has seen two events in the direction 
of the Large Magellanic Cloud (LMC) p[; the OGLE collaboration has seen 12 events in the 
direction of the galactic bulge f^ ; and the MACHO collaboration has seen three events in 
the direction of the LMC and more than 40 in the direction of the galactic bulge |^ . 

The study of microlensing toward the galactic bulge mainly probes the structure of the 
inner galaxy, while microlensing toward the LMC mainly probes the dark halo |jl|, Q. The 
probability that a given star is being microlensed by a foreground object is referred to as the 
optical depth for microlensing (= r). Expectations for the bulge were tbulge — 1 x 10~^, 
largely due to lower-main-sequence stars in the disk 0; OGLE reports an optical depth 
that is about a factor of three larger, tqcle = 3.3 ± 1.2 x 10^^ 0, and the rate observed by 
MACHO may be even higher ||^. Expectations for an all-MACHO halo were tlmc — 5 x 10~^ 
1^ (because data poorly constrain the halo, uncertainties in this estimate are large, almost 
a factor two either way p, ^, |10[)- Based upon 9 million-star years of observations, the three 
events observed, and estimates of their efficiencies (between 20% and 40%) 0], the MACHO 
data indicate that tmacho — 1 x 10~^. The EROS data indicate a similar optical depth [||. 

Evidence that spiral galaxies, including our own, are embedded in extended, dark (roughly 
spherical) halos comes from galactic rotation curves, the study of satellite galaxies and binary 
galaxies, the kinematics of the globular clusters in our galaxy, the warping of galactic disks, 
and the flaring of neutral hydrogen gas associated with disks pTS] , ^, [T^. Galactic halos are 
repositories for both nonbaryonic dark matter (mainly slowly moving particles, or cold dark 
matter, since fast moving particles such as light neutrinos move too fast to be captured) and 
baryonic dark matter. Experimental efforts to detect nonbaryonic dark matter have focussed 
on our own halo |l^. Determining the mass fraction of the halo in baryons is crucial for 



estimating the amount of nonbaryonic matter that may exist in our galaxy. 

Our purpose here is to use the microlensing data to draw conclusions about the MACHO 
fraction of the halo — and from it the fraction of the halo that could be particle dark matter. 
Since the expected microlensing rate in the direction of the LMC depends upon galactic 
modelling 0, ^ |T^ and the LMC microlensing statistics are small, we adopt the following 
strategy. The rotation curve, local projected mass density, distribution of luminous material 
in the disk and bulge, and bulge microlensing rate are used to constrain the halo model. For 
acceptable galactic models we calculate tlmC) and then, based upon the observed LMC rate, 
make inferences about the MACHO fraction of the halo, concluding that it must less than 
about 30%. 



Models of the Milky Way have three major components |T6|, a central bulge, a disk, 
and a spherical halo, with large uncertainties in the parameters that define all three. The 
basic picture of the bulge has evolved from spherical to recent indications that it may be 



closer to a bar |T7|. We follow Dwek et al. [|I^ who have utilized DIRBE surface brightness 



observations to construct a triaxial bulge model: 
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where the bulge mass Mbar = 0.82Mo, the scale lengths a = 1.49 kpc, b = 0.58 kpc and 
c = 0.40 kpc, and the long axis is oriented at an angle of about 10° with respect to the line 
of sight toward the galactic center. The bulge mass is not well determined, and we consider 



Mbar = 1, 2, 3, and 4 x IO^^Mq [y, [Tg]. 



The bulk of the luminous matter in the disk follows a double exponential distribution 
| 20| . There is some evidence that the disk has both a "thick" and a "thin" component |pO|| . 



We take the sum of a "fixed," thin luminous disk, 

PLUM(r, z) = ^ exp[-(r - ro)/r,]e-l^l/\ (2) 

with scale length r^ ~ 3.5 kpc, scale height h = 0.3 kpc, and local projected mass density 
Slum = 25M0pc^^ [21], a "variable" component. For the variable component we use the 



same scale length, but vary the thickness, very thin {h = 0.15 kpc), thin {h = 0.3 kpc) and 
thick {h = 1.5 kpc), and local projected mass density Svar- We also explore a variable 
component whose surface density scales as 1/r rather than exponentially [^. A number 
of studies of stars motions constrain the local projected mass density within a distance 



of 0.3 kpc — 1.1 kpc of the galactic plane [22|. As a conservative bound we require that 



STOT(lkpc) = /_ikpcP(^05 z)dz = 25 — lOOMopc"^; this implies the variable component of 
the disk must satisfy SvAR(lkpc) < 75M0pc~^. 

The third component is the halo. We assume independent isothermal distributions for 
the MACHOs and cold dark matter with core radii Oj = 2, 4, 8, and 16 kpc. 
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where i = MACHO, cold dark matter and po,j is the local mass density of component i. 
More complex halo models are possible; e.g., flattened halos p, 0. We do not expect such 
refinements to significantly alter our basic conclusions; they only serve to increase slightly 
the theoretical uncertainties. 

The average optical depth for microlensing a distant star by a foreground star is P| 

^ 47rG /o°° dspjs) /p dxp{x)x{s - x)/s 
^ & Jo^dspis) ' ^ ' 

where p is the mass density in stars, s is the distance to the star being lensed, and x is 
the distance to the lens EM. In estimating the optical depth toward the bulge, we consider 



lensing of bulge stars by both disk and bulge objects; for the LMC we consider lensing of 
LMC stars by halo and disk objects. We have assumed that the threshold for the detection 
of microlensing is a brightening of 1.34 (which roughly corresponds to the experimental 



thresholds). Further, while the microlensing rate more closely describes what is measured, 
it depends upon detailed knowledge of the velocity distribution of the lenses, and previous 
analyses ||^, |^ have found that optical depth correlates well with lensing rate. 

Kinematic constraints to the galactic model come from the circular rotation speed at our 
position (= Vc) and the requirement that the rotation curve be approximately flat between 
about 4kpc and 18kpc. We adopt the lAU value of 220kms~^ for v^ with an uncertainty 
of ±20kms~^, and we take our distance from the galactic center to be tq = 8.0, 8.5, and 
9.0 kpc. For the flatness constraint we follow our previous work |^ in requiring that the total 
variation in v{r) be less than 14% over the aforementioned range. 

We construct our suite of viable models as follows. Starting with a disk (exponential or 
1/r) with local surface density Svar and a bulge of mass Mbar we compute tbulge, the 
optical depth to Baade's window, galactic coordinates (1°, —4°), and the contributions of 
the disk and bulge to the rotation curve at r = tq. For a choice of halo parameters this then 
determines the local halo density, the full rotation curve, and the optical depth to the LMC 
We deem a model viable if: (a) tbulge > 2.0 x 10~®, (6) the rotation curve is sufficiently 
fiat, and (c) tlmc is in the range 0.5 — 2.0 x 10^^. The last condition primarily constrains 
the baryonic mass fraction and does not eliminate many models. 

Our results are summarized in Figs. 1-4. We find that in general the disk alone does 
not provide sufficient lensing to explain the event rate seen toward the bulge. While tbulge 
increases with Svar, ^var reaches its upper bound of 75Mq pc~^ or the rotation-curve 
constraint is violated before tbulge = 2 x 10^®. For fixed Svar the disk density along the 
line of sight is about the same, and the difference between a thick and thin disk is negligible. 
However, there is less mass in a thin disk and a larger Svar is permitted without violating 
the rot at ion- curve constraint. A very thin disk [h = 0.15 kpc) does not help further: the 
line of sight to Baade's window passes above most of the disk material. The results for disks 
with a 1/r density distribution are not significantly different from those of exponential disks. 



The bar is a much more efficient source of lensing [^ . In all viable models a bar mass of 
at least 2 x lO^^M© is required; if tbulge is determined to be greater than 3 x 10~^, as may 
well be the case, a bar mass of at least 3 x IO^^Mq seems unavoidable. Should the optical 
depth be 4 X 10~^, a bar mass of 4 x IO^^Mq may be indicated. We note that even higher 
bar masses (Mbar > 5 x IO^^Mq) make it difficult to achieve a fiat rotation curve interior 
to the solar radius and thus are not viable. 

Turning now to the optical depth to the LMC (Fig. 2), we find as expected that a thin 
disk makes a negligible contribution to tlmc? while a thick disk can provide a significant 
contribution (toward the LMC r oa h). In fact, a model with a heavy bar and a thick disk 
can account for both the bulge and LMC rates without recourse to MACHOs in the halo. 
We also mention that the LMC microlensing rate is small enough that a significant part of 
it (~ 0.5 X 10~^) could be due to microlensing by objects in the LMC itself ||2^ . 



Our most striking result is that all viable models have a significant halo component. This 
can be traced to the fiat-rotation curve constraint (even though it is very conservative) and 
the upper limit to local projected mass density |]23|, and is essentially independent of the 
bulge optical depth. That is, while tbulge provides us with information about the disk and 



bulge components of our galaxy, the halo parameters are relatively insensitive to it. This can 
be seen in Fig. 3, where the local halo density in acceptable models is around 5 x 10~^^ g cm~^ 
with more or less the same uncertainty as previous estimates (see e.g., Refs. p7[|). 

Even more than the uncertainty in the observed LMC microlensing rate, imprecise knowl- 
edge of the galactic model dominates the uncertainties in determining the MACHO fraction 
of the halo. However, in most viable models, an all- MACHO halo results in an optical depth 
to the LMC that is many times that observed (Fig. 2), and a significant nonbaryonic halo 
component seems indicated (Figs. 3 and 4). The viable models with the largest MACHO 
fraction have thin disks, small bar masses, large core radii and high local circular velocity. 

To summarize: (1) A single component of the galaxy cannot account for both the bulge 
and LMC events. For example, a spherical halo predicts an LMC rate that is slightly higher 
than the bulge rate; a thin disk cannot account for either the bulge rate or the LMC rate. 
A thick disk can explain the LMC events, but not the bulge events. (2) The most promising 
model for explaining the high bulge rate is an asymmetric bulge (bar) that lenses itself 



25| with a lesser contribution from a thin disk. (3) Viable models of the galaxy have a 
significant halo component and the LMC rate expected for an all-MACHO halo is many 
times that observed. (4) While the present data cannot preclude an all-MACHO halo, it 
appears that the fraction of the halo in MACHOs is less than 30%. This is consistent with 
searches for faint halo stars that indicate that their contribution to the halo mass is small 



(provided that the mass function of halo stars is "smooth;" see e.g., Ref. |^ 



If the bulk of the halo is not in the form of MACHOs, what is it? While it is not 



impossible that it could be baryonic, in a more diffuse form, e.g., clouds of neutral gas [^ 
cold dark matter is a more compelling possibility. In an fi = 1 cold dark matter model, the 
naive expectation for the baryonic mass fraction in our galactic halo is Jb = ^b- Based upon 
primordial nucleosynthesis 0.01/t~^ < f^B < 0.02/i~^ [^, which implies Jb — 0.04 — 0.2. The 
expected baryon mass fraction for models with an admixture of massive neutrinos [Q^, ~ 
0.2 — 0.3) is only slightly higher. For cold dark matter models with a cosmological constant, 
Q\ ~ 0.8, fs — 0.1 — 0.2. If the baryonic halo has undergone moderate dissipation, the 
baryonic mass the baryon fraction in the inner part of the galaxy can be increased, though it 
is still expected to compose less than half of the local dark matter density [^]. From Fig. 4, 
it is clear that the baryonic mass fraction in our halo implied from microlensing is consistent 
with any of these scenarios which provides further motivation for the ongoing experimental 
efforts to directly detect neutralinos and axions in our own halo. 
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Figure Captions 

Figure 1: Optical depth to the bulge for bar masses of 1, 2, and SxIO^^Mq (bottom to top) 
and thick (broken) and thin (solid) disks as a function of the local projected mass density 
SvAR for models that satisfy the kinematic constraints (ro = 8.5 kpc and Vc = 220kms~^). 

Figure 2: Optical depth to the LMC from an all- MACHO halo (upper lines) and thick 
(broken) and thin disks (solid) for bar masses of 1, 2, and SxIO^^Mq (right to left) as 
a function of Svar for models that satisfy the kinematic constraints (rg = 8.5 kpc and 
Wc = 220kms"^). 

Figure 3: Distribution of local cold dark matter density in viable models for Svar = 5 — 15, 
25-35, 45-55, and 65 — 75M0pc~^. Since the halo MACHO fraction in most viable models 
is small, the local cold dark matter density is approximately equal to the total local halo 
density. 

Figure 4: Distribution of halo MACHO mass fraction in viable models for Svar = 5 — 15, 
25-35, 45-55, and 65 - 75M0pc-2. 
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